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Abstract 

We develop a new approach for the Anderson localization problem. The implementation of 
this method yields strong numerical evidence leading to a (surprising to many) conjecture: The 
two dimensional discrete random Schrodinger operator with small disorder allows states that are 
dynamically delocalized with positive probability. This approach is based on a recent result by 
Abakumov-Liaw-Poltoratski which is rooted in the study of spectral behavior under rank-one 
perturbations, and states that every non-zero vector is almost surely cyclic for the singular part 
of the operator. 

The numerical work presented is rather simplistic compared to other numerical approaches in 
the field. Further, this method eliminates effects due to boundary conditions. 

While we carried out the numerical experiment almost exclusively in the case of the two 
dimensional discrete random Schrodinger operator, we include the setup for the general class of 
Anderson models called Anderson-type Hamiltonians. 

We track the location of the energy when a wave packet initially located at the origin is evolved 
according to the discrete random Schrodinger operator. 

This method does not provide new insight on the energy regimes for which diffusion occurs. 
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1 Introduction 

In 1958 P. W. Anderson [Anderson 1958 suggested that sufficiently large impurities in a semi- 
conductor could lead to spatial localization of electrons, called Anderson localization. Although many 
physicists consider the problem solved, many mathematical questions with striking physical relevance 
remain open. The field has grown into a rich mathematical theory (see [Germinet-Klein-Schenker 2007 
and |Ghribi-Hislop-Klopp 2007] for the study of different Anderson models; for refined notions of An- 



derson localization see del Rio-Jitomirskaya-Last-Simon 1986| and [Last 2007] ). 



We consider the discrete random Schrodinger operator in dimension d which is given by the self- 
adjoint operator 

H u = - A + < • > s i > 5 i on l2 ( zd )- 

Here Si G / 2 (Z d ) assumes the value 1 in the i— th entry, i = (ii, Z2, . . . , id), and zero in all other entries 
(see equation Q below for an example in two dimensions). The random variables uji are i.i.d. with 
uniform distribution in [— c, c], i.e. according to the probability distribution P = (2c) -1 Hx r _ A x - 

[ c,cj 

The Laplacian describes a crystal with atoms located at the integer lattice points Z d . Adding the 
random part can be interpreted as having the atoms being not perfectly on the lattice points, but 
randomly displaced. 

This paper pertains to one of the "weaker" definitions of localization which is equivalent to the 
property that the spectrum of the operator is almost surely purely singular. 

It is known that for d = 1, Anderson localization is produced by random disorders of any strength 
and at all energies (analytic proof, see e.g. |Cycon-Froese-Kirsh-Simon 1987 , |Carmon~ Lacroix 1990] . 
or |Figotin-Pastur 1991| ). 

For d > 2, localization is proved analytically for disorders c above a certain dimension-dependent 

threshold q. The first result of this type can be found in IFrohlich-Spencer 1983 . Simpler proofs 
-i -i , , , -i r i • rr~' * t i i '-.^oi r; 1 . ' -.^.t^ i 



and better constants can be found in |Aizenman-Molchanov 1993] as well as [Simon 1994B] , 

Diffusion may hence only occur for small disorder c. This is precisely the question we are addressing: 
Does the two dimensional discrete random Schrodinger operator exhibit Anderson localization for small 
disorder with non-zero probability? 



*The author is partially supported by the NSF grant DMS-1 101477. 
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The main contribution of this paper is the introduction of a new numerical approach, and its 
implementation for the two dimensional discrete random Schrodinger operator. This application 
supports the following conjecture which overthrows the widely spread belief that localization takes 
place for random disorders of any strength. 

Conjecture 1.1 (Derealization conjecture). For disorder c < 0.7, the two dimensional discrete 
random Schrodinger operator does not exhibit Anderson localization with positive probability; in the 
sense that it has non-zero absolutely continuous spectrum with positive probability. In particular, we 
do not have dynamical localization with positive probability for small disorder. 

While this conjecture is based on deep analytical results, we included an explicit statement of 
the main tool, see Corollary |3.3[ so that the applied numerical sections [4] through [7] are essentially 
self-contained. We would like to point out the important feature of this Corollary which makes a 
numerical experiment feasible: It suffices to track the evolution (under the random Hamiltonian) of 
just one vector! 

In Section [2j we introduce Anderson- type Hamiltonians, a very general notion of Anderson model 
which includes many of those studied in literature. The methods described within this paper can 
be extended to most Anderson- type Hamiltonians. We further explain the notion of singular and 
absolutely continuous spectrum, as well as the corresponding parts of the operator. In Section [3] 



we state an improvement, Theorem 3.1, of a result by Jaksic and Last concerning the cyclicity of 
vectors for the Anderson-type Hamiltonian. We further explain how analytic results are used to prove 
the main tool, Corollary |3.3| Section [4] is devoted to a description of the numerical experiment. A 
summary of numerical results and the conclusions can be found in Section [5] In Section [6] we verify 
the performance of the method in many examples, e.g. for large disorder, for the free/unperturbed 
two dimensional discrete Schrodinger operator, for the one dimensional discrete random Schrodinger 
operator. We further include an investigation of the distribution of energies after repeated application 
of the random operator of a wave packet initially located at the origin. We briefly remark on computing 
and memory requirements in Section [7[ 

Acknowledgements. The author would like to thank A. Poltoratski for suggesting the initial 
mathematical idea of the experiment, and G. Berkolaiko for the insightful discussions concerning many 
aspects of this research as well as for reading and making useful comments on most parts of this paper. 
Further, she would like to thank J. Kuehl for running initial experiments using a code of his, as well 
as for being such a wonderful husband. 



2 Preliminaries 

2.1 Anderson-type Hamiltonians, discrete Schrodinger operator 

While the numerical experiment within pertains to the discrete random Schrodinger operator, we 
define so-called Anderson- type Hamiltonians which were first introduced in [Jaksic ^Last 2000] . The 
advantage of making this general definition is that this notion is a generalization of many Anderson 
models discussed in literature. In particular, the method described within this paper can be applied 
to many other interesting Anderson models. 

For n G N consider the probability space Vt n = (R, S, /i n ), where B is the Borel sigma- algebra on 
R and \i n is a Borel probability measure. Let Vt = n^Lo be a P r °duct space with the probability 
measure P on Q introduced as the product measure of the corresponding measures on Q n on the 
product sigma-algebra A. The elements of ft are points in R°°, uo = (lji,lj2, - .-) £oy u n e Q n . 

Let H be a separable Hilbert space and let {(p n }neN be a countable collection of unit vectors in %. 
For each uo G ft define an Anderson-type Hamiltonian on T~L as a self-adjoint operator formally given 
by 

Hco = H + Vu,, where V u = ^ oj n < • , (p n > <p n . (1) 
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Except for degenerate cases, the perturbation V u is almost surely a non-compact operator. It is 
hence not possible to apply results from classical perturbation theory to study the spectra of if^, see 
e.g. |Birman^Solomj ak 1986 and |Kato 1980 . 



In the case of an orthogonal sequence {<^ n }, this operator was studied in [Jak sic-Last 2000| and 
[Jaksic-Last 2006^ . Probably the most important special case of an Anderson-type Hamiltonian is the 
discrete random Schrodinger operator on l 2 (Z d ) 

Hf(x) = - A f{x) = J2 (/(* + U ) ~ fW)> With = = { J X el Z Ue ( 2 ) 

M=i 1 



2.2 Singular and absolutely continuous parts of normal operators 

Recall that an operator in a separable Hilbert space is called normal if T*T = TT*. By the spectral 
theorem operator T is unitarily equivalent to M z , multiplication by the independent variable z, in a 
direct sum of Hilbert spaces 

H = ® J H(z) dfi(z) 

where \i is a scalar positive measure on C. The measure \i is called a scalar spectral measure of T. 

If T is a unitary or self-adjoint operator, its spectral measure \i is supported on the unit circle or 
on the real line, respectively. Via Radon decomposition, \i can be decomposed into a singular and 
absolutely continuous parts /i = /i s + /i a c- The singular component /i s can be further split into singular 
continuous and pure point parts. For unitary or self-adjoint T we denote by T ac the restriction of T to 
its absolutely continuous part, i.e. T ac is unitarily equivalent to M t \^ ju^d^ (t)' Similarly, define the 
singular, singular continuous and the pure point parts of T, denoted by T s , T sc and T pp , respectively. 



3 Theoretical background and the main tool 

Let us explain how the main theoretical tool used to indicate derealization, Corollary |3.3[ is deduced 
from the following theorem. Let us remind the reader that we use the term derealization to mean 
the existence of absolutely continuous spectrum almost surely, and that such derealization implies 
dynamical derealization. 

A sequence {ip n } C T~L is called a representing system, if every vector ip G T~L can be represented as 
a series ip = ^a n (p n that converges with respect to the norm of %. Note that, bases are representing 
systems. However, unlike in the case of a basis, a representation of a vector need not be unique. 

We use the following Theorem, see [Aba kumov-Liaw-Poltoratskii 2012] . 

Theorem 3.1. Let H u be the Anderson-type Hamiltonian introduced in equation 0. Suppose that 
the probability measure P is a product of absolutely continuous measures and {(p n } is a representing 
system in H. Assume that there exists a vector ip G T~L that is cyclic for H^, P- almost surely. Then 
any non-zero cp G T~L is cyclic for H^, P -almost surely. 

It is well-known that if an Anderson-type Hamiltonian is purely singular almost surely then 
it is cyclic almost surely. Equivalently, if such an operator is not cyclic with positive probability, 
then there are energies which are diffusive with non-zero probability. A proof of almost-sure cyclic- 
ity of the singular part (H^)* and almost-sure cyclicity of certain specific vectors can be found in 
[Jaks ic-Last 2QQ6j and for the discrete Schrodinger operator in [Simon 199 4A j. Together with the 
latter theorem and the fact that the S n form a basis of l 2 (Z d ) we obtain the following result. 

The following two statements were formulated by A. Poltoratski (private communications). 

Corollary 3.2. Assume the hypotheses of Theorem \ 3.1\ If Anderson localization occurs, then the 
orbit of any non-zero vector under the operator H u is almost surely dense in the Hilbert space T~L. 
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In other words, fix ^ / G T~L. If has purely singular spectrum almost surely, then for all 
v G T~L with norm \\v\\^ = 1 the distance from the vector v to the span of the orbit of f under the 
operator is zero with respect to F-almost surely, i.e. with respect to F-almost surely 

dist(>, span{^/ : k G {0, 1,2,..., n}}) 0. 

More specifically, we will apply the following immediate consequence. 

Corollary 3.3. Consider the discrete random Schrodinger operator given by equations and ^ 
in dimension d = 2. Let Ui, i G Z 2 ; be i.i.d. random variables with uniform (Lebesgue) distribution 
on [— c, c], c > 0. To prove derealization (i.e. the existence of absolutely continuous spectrum with 
positive probability), it suffices to find c > for which the distance 

dist(£n, spanji^Joo : k = 0, 1, 2, . . . , n}) as n — > oo (3) 

with positive probability. 

Remark 3.4. The converse of Corollary \3.3\ is not true. Hence we cannot draw any conclusions, if 
the distance between a fixed (unit) vector and the subspace generated by the orbit of another vector 
tends to zero. In particular, we cannot conclude that there must be localization. Even if we show (|3| 
for many or c alV vectors (instead of just Sn), it could be possible that the absolutely continuous part 
has multiplicity one and that 5qq is cyclic, that is, / 2 (Z 2 ) = clos spanji^^oo : k G N U {0}}. 



4 Method of numerical experiment 

Let us explain the computational approach used to indicate diffusion. 

Consider the discrete Schrodinger ope rato r given by and ([2| with random variable uo distributed 
according to the hypotheses of Corollary 3.3 Fix the vectors 8qq £ ^ 2 (Z 2 ) and 5u G / 2 (Z 2 ), that is 
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Notice that 



DZ c := dist^u, spanfff^oo 



0,1,2, 



simply describes the distance between the unit vector Sn and the subspace obtained taking the closure 
of the span of the vectors <S o> #uAo, H%5 o, . . . , H™5 00 . 

In virtue of Corollary 3.3 , we obtain derealization, if we can find c > for which Q happens 
with non-zero probability. 

In the numerical experiment, we initially fix c and fix one computer- generated realization of the 
random variable uo (with distribution in accordance to the hypotheses of Corollary 3.3). We then 
calculate the distances c for n G {0, 1, 2, . . .}. In Subsection 
approach used to compute D™ . 



4.2 



(below) we describe the numerical 



Assuming that we know D™ c for n = 0, . . . , 4500, let us find a lower estimate for the limit 



D U)C := lim D" . 
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0.8875F 




0.8865 



Figure 1: Typical trend for the distance D™ c as a function of n. Notice the fine vertical scale. 



Figure [T] shows typical trends for the distance D™ c as a function of n. As the first n = 200 points 
do not contribute to the value of this limit and were generally rather irregular, we have omitted those 
points in what follows. (On the side, notice that the graph shows a decreasing function, as expected.) 

While those results at this point looked fairly promising, they were not yet satisfactory. Most of 
all they do not provide a reliable estimate for the limit D^ c . In order to obtain such an estimate for 
A^,c? we re-scaled the horizontal axis in Figure [l] by a negative power n~ a (power of the reciprocal, 
so that the horizontal axis is reverted) and approximated the resulting graph by a line. 

The re-scaled graph is shown in Figure [2] Subsection |4.3| contains information about the choice of 
the re-scaling factor and explains why, for appropriately small disorders, the graph does not decay to 
zero, e.g. logarithmically. The subtleties of choosing the re-scaling parameter a are the reason why 
we do not expect derealization with probability one in the Derealization Conjecture but rather 
with non-zero probability. This decision is explained further in Subsection |4.3| 

The value of D u c is estimated by the intercept y UiC of the approximating line. 

Since the re-scaled graphs in Figure [2] were sometimes rather noisy (e.g. a line through the steeper 
sections of the graph has a lower y— intercept), we decided to include a lower estimate L UyC for y^^ 
given by the minimum intercept of the lines passing through any two consecutive points. 

Summarizing the last few steps, we have 



(5) 



Finally, we repeat the experiment for many values of c and many computer-generated realizations 
for the random variable uj. Concerning the different realizations, throughout we took the minimum of 
y UyC and L UyC over all the different computer-generated realizations of uj. Roughly the goal is to show 
that for some c > 0, the limits D W)C are bounded away from zero for many realizations uj. 
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n- 2 

Figure 2: Re-scaling of the horizontal axis in Figure [I] using the best exponent a = 0.2. The 
y— intercept of the approximating line is the estimate y UyC of the value for D W;C . 

4.1 When to fix the realization uj 

In the experiments described here, we had fixed c and uo at the beginning. For fixed c = 0.1, we 
also computed several cases for which we chose a different realization uo each time we applied he 
random operator. In other words, for c = 0.1 we fix countably many realizations of uo each indepen- 
dently distributed and each in accordance with Corollary |3.3| Let those realizations be denoted by 
uji, i G N. Then we compute the distance between Sn and the closure of the span of the vectors 
^oo, H OJl Soo, H OJ2 (H UJl Soo), H UJ3 (H OJ2 H UJl Soo), etc. 

The results obtained from this setup agreed very well with the ones described in Section [5] below. 

4.2 Computing the distance D™ c 

For fixed n and uo let us briefly explain the computational approach to obtain D™ c . The main idea is 
to apply the Gram-Schmidt orthogonalization process in order to recursively compute c . 
Take mo = #oo and D® = 1 (since #oo and Su are orthonormal) . 

In order to compute D™^ 1 , assume we have an orthonormal basis {mo,mi,m2, . . . ,m n } for the 
linear subspace 

X n := span{#^oo : k = 0, 1, 2, . . . , n} of I 2 (I?). 

Let us find an orthonormal basis for X n +\. According to the Gram-Schmidt orthogonalization 
process, we define m n+ i to be the unit vector in the direction of 

n 

Hytrin - ^ < H^rrin, mi > m t . 

1=0 

The following proposition says that all but the two last terms in the sum are zero. We learned this fact 
and its proof from a conversation with M. Hastings. This simplification reduces the required memory 
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by the order n (from 0(n 3 ) to 0{n 2 )). 

Proposition 4.1. The vector H^rrin is orthogonal to mi for all I = 0, 1, 2, . . . , n — 2. 

Although this result seems to be well-known to the physics community, we include the short proof 
by mathematical induction on n. 

Proof. Consider n = 2. Assume that we have computed the orthonormal vectors mo, mi and m^ 
(via the Gram-Schmidt orthogonalization process). Since the operator is self-adjoint, we have < 
H UJ m2^ 171 o >=< m2, H^nio >. Since mi = H^rriQ— < H^mo^mo > mo and because m2 is orthogonal 
to mo and m<i, we obtain 

< i^ u; m2,m >= (m 2 ,mi+ < H^rriQ, m > m ) = 0. 

Assume that the statement of the proposition is true for some n — 1 > 2. It remains to show that 
the statement is true for n. Assume that we have computed an orthonormal sequence mo, mi, . . . , m n . 
For I < n — 2, it suffices to show that < H u m ni mi >= 0. By following the argument for the base case, 
we obtain < H u m n ,mi >=< m n , H^rrij >= (m n ,m^ + i+ < H u mi, mj > m{). The latter expression 
equals zero, because we assumed I < n — 2 and the orthogonality assumption on mo, mi, ... , m n . □ 

According to the latter proposition we take m n+ i to be the unit vector in the direction of 

m n+ i mn+1 ,, , where m n+i H^m n - < H u m n ,m n -i > m n _i- < R^m n ,m n > m n . 

Now, the distance D^ 1 of the vector Sn to the subspace X n+ i equals the Euclidean norm 

= K+1H2 with e n+ i = 5 U - Pn+iSu, (6) 

and where P n +i denotes the orthogonal projection from / 2 (Z 2 ) onto X n +\. 

A little more analysis allows us to simplify the latter expression. The following expression is closely 
related to the dimensionless scaling parameter that occurs in the so-called Thouless criterion. 

Proposition 4.2. We have (D™^ 1 ) 2 = 1 — Y^i=2^ x i) 2 w ^ere xi denotes the (1, 1) — entry of mi. 

Proof. By the definition of the m n 's, the vectors are recursively given by 

e = £11 and e n+i = e n - < e n ,m n+ i > m n+i . 

We use this same recursive definition for e n in the inner product and the fact that the m n form 
an orthonormal sequence to obtain 

e n +i e n - < (e n _i- < e n _i,m n > m n ),m n+ i > m n+i = e n - < e n _i,m n+ i > m n+ i. 

Repeated application of this argument yields that 

e n +i = e n - < e ,m n+ i > m n+ i e n - < 5u,m n +i > m n+i . 

We proceed to replace e n by its recursive definition and so on, until we obtain 

n+l 

e n +i = *n - ^2 < ^ n ' m ^ > mi ' 
1=2 

The proposition follows from equation (|6| and the Pythagorean Theorem, since the mi form an 
orthonormal sequence and all of them are orthogonal to e n +i. Also notice that ||#ii||2 = 1- □ 
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4.3 Choice of the re-scaling parameter 

For each fixed c and cj, the re-scaling exponent a is chosen so that the re-scaled graph of the distance 
function (see Figure [2| satisfies the least square property; that is, the error when approximating the 
graph by a line is minimal. With this exponent we then find the corresponding linear approximation 
for the re-scaled distance function. 

We include an extract of the table of best re-scaling exponents a which satisfy the least square 
property for our data. As the values for y u ^ c were not very sensitively dependent on the precise value 
of a, we used a rather coarse mesh a = 0.05 : 0.05 : 0.85 and refined using a = 0.01 : 0.01 : 0.05, if the 
best re-scaling exponent was below 0.05. Each entry in the table corresponds to a different realization 
of the random variable uj. 



c 


.1 


.15 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


.9 


1 


1.2 


1.3 




0.2 


0.04 


0.2 


0.25 


0.1 


0.15 


0.1 


0.1 


0.05 


0.02 


0.05 


0.1 


N/A 




0.1 


0.2 


0.25 


0.1 


0.04 


0.1 


0.1 


0.05 


0.05 


0.02 


N/A 


N/A 


0.05 




0.2 


N/A 


0.2 


0.2 


0.15 


0.1 


0.2 


0.1 


0.1 


0.04 


0.03 


N/A 


N/A 




0.05 


0.1 


0.05 


0.15 


0.15 


0.1 


0.05 


0.1 


0.05 


0.03 


0.05 


N/A 


N/A 



(7) 



The N/A indicates that for this particular realization, even the re-scaling parameter a = 0.02 
yields a concave graph. For this realization, we do not obtain any information. For values of c > 1.2 
many realizations did not yield a reasonable best fit parameter a. No statement can be made for such 
disorders. 

Since we can only investigate finitely many randomizations, and one of the realizations for a fairly 
small value of c = 0.15 yielded an inconclusive result, we decided to conjecture derealization with 
non-zero probability in the Derealization Conjecture |1.1[ rather than almost surely. 

The existence of a positive re-scaling factor implies that the graph in Figure [I] will not decay to 
zero, e.g. logarithmically. Indeed, if we use a re-scaling factor smaller than the one in the table will 
result in a 'globally concave' graph for the distances D™ c . In this case, the y— intercept of the line 
lies below the value expected for D^ c . 

5 Conclusions 

As mentioned in Section |4j for a fixed c we chose many realizations uj. We took the minimum of the 
resulting quantities for L^^ and 2/ WjC (the intercept of the approximating line and the minimum 
y— intercept of the lines passing through any two consecutive points, respectively). 

Figure [3] shows L UyC and y UiC as a function of c. Being rather cautious, we say that a negative 
value for y UyC indicates that the orbit of Soo may not span the whole space. Hence the final conclusion 
of this numerical experiment is precisely the Derealization Conjecture |1.1| 

As it was explained in the remark following Corollary |3.3[ we cannot conclude localization even if 
Voj,c < 0- Therefore, the experiments do not imply localization for larger values of c. 

6 Further supporting the credibility of the method and the 
numerical experiments 

Apart from the usual tests (the program is running stably, checking all subroutines, many verifications 
for small n) , we have also tested the code and versions for other models: the free/unperturbed two 
dimensional Schrodinger operator and the one dimensional random Schrodinger operator. We briefly 
summarize the results, in order to provide verification for the correctness of method and code. 

Further, we provide information of the energy distribution in terms of the distance from the origin 
of the evolution of the vector $00, describing how a wave packet which was initially located at the 
origin changes as time progresses. 
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Figure 3: As a function of c we show y u ^ c (Vs; larger values) and L u ^ c (circles; smaller values). Notice 
that D^^c ~ Voj,c > L^^c > for c < 0.7, supporting our Delocalization Conjecture 



6.1 Free discrete two dimensional Schrodinger operator 

When we apply the free discrete Schrodinger operator H = Ho to the vector #oo, it immediately 
becomes clear that Hdoo as well as all vectors H n 5oo, n G N U {0}, are symmetric with respect to the 
origin. In dimension d = 2, it is not hard to see that the distance between Sn and the orbit of ^oo 
under H is at least a/3/2 w 0.8660. Indeed, we have 

dist(Jn, closspan{iTM o : n G NU {0}}) > mmdist(u x ,5 n ) = V3/2, 



where 
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In the experiments for the free discrete two dimensional Schrodinger operator we obtained a 
y— intercept of the approximating line approximately equals 0.8867. The re-scaled graph of distances 
still had a very convex shape, so the actual distance as n — » oo would be bigger. In fact, we have 
extracted from Figure [4] an upper estimate of 0.8868 by zooming in. Therefore, the distance must lie 
in the interval [0.8866,0.8868]. 
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0.8868 




Figure 4: Convex shape for the distance from Su to the orbit of Joo under the free Schrodinger 
operator in two dimensions. The approximating line has y— intercept y UyC « 0.8866. 



6.2 Verifying localization for the one dimensional random Schrodinger op- 
erator 

Consider the discrete random Schrodinger operator in one dimension, see e.g. equations and Q 
with d = 1. For this operator, it is well known that localization occurs for random disorders of all 
strengths (in particular, for small values of c) and at all energies. 

We have adopted and applied this computational approach for the discrete random Schrodinger 
operator in one dimension. Figure [5] shows a typical re-scaled graph of the distance 

Dl c = dist^i, span{i^£ : k = 0, 1, 2, . . . , n}) for n = 3000, 3001, 3002, 15000 

for the disorder c = 0.05. With a re-scaling exponent of a = 0.09, the graph of D™ c is still concave, 
so that the y— intercept of the approximating line is an upper estimate of D W;C the limit. Therefore 
we have 

D u , c < y u ,c = -0.0543. 

While we know by the remark following Corollary |3.3| that this experiment does not allow us to 
conclude that there is localization, the result still provides support for the credibility of the method 
at hand as well as the numerical design. 



6.3 Diffusion of energy for small values of c 

We present the distribution of energies of a wave packet initially located at the origin as the random 
operator is repeatedly applied. By distribution of energies, we mean how much of the energy is located 
at which 'distance' from the origin. 

For example, in order to obtain how much energy of the vector rrtk (defined in Subsection 4.2) is 
at 'distance' 2 from the origin, we use the elements of which are located on the diamond for which 
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°' 6 (^42 0.43 0.44 0.45 0.46 0.47 0.48 0.49 

-.09 



Figure 5: Discrete random Schrodinger operator in one dimension with disorder c - 
we have y U}C « —0.0543 and a convex graph. 



0.05. For a = 0.09, 



the matrix 




has entries equal to 2. The energy E(2, k) of the vector at 'distance' 2 from the origin is equal to 
the Euclidean norm over the elements in this diamond. In general, we have 



E(l,k) 




(8) 



for the energy E(l,k) of the vector m\~ at 'distance' I from the origin. Here (rrik)ij refers to the 
entry of the 2 x 2— matrix m&. 

By small modifications of our programs, we have extracted the location of the energy the vector 
Sqq evolves under the random Hamiltonian for the values c = .1, c = 1 and c — 5 of disorder, see 
Figure [5J In accordance with our Derealization Conjecture the energy for small disorder is far 
away from the origin whereas it is concentrated close to the origin for large disorder. 

Figures [6] and [7] shows the energy distribution of H^Soq for n = 2999 for values of c ranging from 
c = .1 to c = 1. Again, the fact that the energy for small disorder is far away from the origin whereas 
it shifts much closer to the origin as the disorder increases, supports Derealization Conjecture 

Both figures are the averages obtained from two realizations for each value of c. And again, 
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diamond n 



Figure 6: For c = 0.05 we show E(l, n), i.e. the evolution of the energy distribution of H™ c #oo fc> r the 
diamonds at distance I from the origin. Notice that the energy travels far out from the origin (the 
diagonal is the farthest possible). 




diamond 



Figure 7: The analog of Figure [6] for c = 0.6 (but with different color scale). The energy remains 
much closer to the origin. 
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0.2 

distance from origin 



Figure 8: The figure shows E(l, 2999) (the Euclidean norm on the diamonds of 777,2999 given by equation 
as a function of the 'distance' I from the origin for several values of disorder c. 



many repetitions of these experiments for smaller values of n were carried out, and the figures shown 
represent the behavior obtained in all repetitions. 

6.4 Precision 

The results are not a phenomenon of numerical errors (e.g. round off errors that sum up over time). 
Indeed, we compared our results with those of a double precision computation. The results agreed 
very well. 

7 On computing and memory requirements 

The implementation uses memory rather efficiently, so that the numerical experiments were mainly 
limited by the length of the computation. On the rather small machines available to us, it took 8 1/2 
hours to complete one realization for one value of c. Since we need to include many realizations of the 
random variable and many values of c, it took even several units a considerable time to finish all the 
computations. 

In order to compute D™ described in Section [IJ our code requires order n 2 (i.e. 0(n 2 )) memory. 
Indeed, in order to carry out the Gram-Schmidt orthogonalization process described in Subsection 
4.2 , we must store matrices of size 0{n) x 0{n). The corresponding code for the d— dimensional 



discrete random Schrodinger operator will require memory size of order 0(n ). 

The random Schrodinger operator on the (dyadic) tree uses memory of 0(2 n ). With the resources 
available to us, memory restrictions would only allow us to compute up to n « 27 for the tree. In 
this case, we cannot produce sufficient data to support the fact that the discrete random Schrodinger 
operator on the tree does indeed exhibit derealization. 
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